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Abstract. - 

We theoretically investigate the mechanical effect of the light-induced dipole-dipole interac- 
tion potential on the àtoms in a Bose-Einstein condensate. We present numerical calculations 
on the magnitude and shape of the induced potentials for different experimentally accessible 
geometries. It is shown that the mechanical effect can be distinguished from the effect of 
incoherent scattering for an experimentally feasible setting. 



Introduction. - The interaction of a radiation field with an ensemble of àtoms has been 
investigated in many distinct contexts since the early works of Lorentz and Lorenz [1] of light 
in dense media. Within one century, researchers came up with an innumerable variety of 
light-matter interaction types [2] and reached a recent summit with the development of làser 
cooling [3] . The dipole character of an atom, driven by an electromagnètic wave, can be used 
to apply forces on àtoms by intensity gradients of light fields [4] but it can also produce forces 
between àtoms [5] . The implication of interacting àtoms on the refractive index was treated 
theoretically for dilute Bose gases [6] with an extended Lotentz-Lorenz model. The effect of 
light-induced or excited state collisions shows up in làser cooling by losses in magneto-optical 
traps [7]. The coherent part of this interaction is recently discussed in basically two directions. 
The one focuses on the properties of the light absorbed and emitted in dense media, namely 
radiation trapping [8] , level shifts [9] , dipole blockade effects as a tool for quantum information 
processing [10] and collective effects such as superradiance [11]. The other main focus tends 
towards new forms of interaction in Bose-Einstein condensates as effective 1/r potentials [12] 
or as a generator for rotons [13]. Recently, the effect of the resonant dipole-dipole interaction 
in a cold cloud of highly excited Rydberg states was observed by broadening of spectroscopic 
lines [14] and by resonant energy transfer between different Rydberg states [15]. 

In this letter, we investigate new physical aspects that arise in dense cold atòmic samples 
irradiated by a near resonant làser beam. Àtoms exposed to an electromagnètic wave responds 
as damped harmònic oscillators and exhibit an alternating elèctric dipole moment. The inter- 
action energy of such dipoles can exceed the one of magnètic dipoles in atòmic ground states 
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by several orders of magnitude. We propose and theoretically investigate an experiment which 
allows to study the coherent interaction of làser induced elèctric dipoles by transferring the 
interaction energy among the dipoles into kinetic energy, which can be probed with Standard 
time of flight techniques. Initially, the dipoles are generated for a certain flash time by a 
làser beam with linear polarization in a spin polarized sample of cold àtoms. During the flash 
time, the dipole moments reach a steady state and the light-induced dipole-dipole interaction 
potential is build up. As a first step, we calculate this potentials for a certain density distri- 
bution including the retardation effects of the dipolar fields and the driving electromagnètic 
wave. The flash time is chosen long enough that the àtoms can evolve in the induced potential, 
namely to gain momentum, but short enough not to change the initial density distribution. 
As a next step we discuss the increase of the initial momentum distribution for different ge- 
ometries of the atòmic cloud and as a function of the angle of the linear polarization of the 
làser light. Finally, the results are compared with parasite effects like spontaneous scattering. 

Outline of the calculation. In the following, we are dealing with oscillating dipoles 
driven by an electromagnètic wave with linear polarization and wave vector k. We assume 
that all dipole moments are of equal strength, oriented parallel and oscillate with the same 
frequeney. The phase between two dipoles depends on the position of the àtoms with respect to 
the electromagnètic wave phase fronts and the interatomic distance. The retarded interaction 
potential for two interacting dipoles with one dipole located at the origin and the other at ro 
reads [16] 
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where d is the absolute value of the dipole moment and fc the wave vector of the driving 
field. Finally, we want to look at a system of N pairwise interacting dipoles with a density 
distribution n(r). The superposed potential for a dipole at position fó is given by 

V dd(?o) = I V AA {f 0l r)n{r)d\ . (2) 

Rcplacing V AA (fo, r) by d 2 cos(/c • (fó — r))y' dd (fó — r), it is possible to rewrite the integral 
as a convolution 

^dd(^o) = / d 2 cos(fc ■ (ro - r))V^ d (ro - ?)n{r)d\ (3) 

to which the convolution theorem can be applied. Since there exists no analytical solution, 
we discreticize the integral on a simple cúbic lattice and use FFT algorithms for evaluation. 

The induced dipole potential V^y changes the momentum distribution of the atòmic cloud. 
In the following, we assume a density distribution of a Bose-Einstein condensate with all 
àtoms having the same phase. The undisturbed wave function can be written as ipo(f) = 
e^ (í V"o(r). The wave function is an eigen-state of the unperturbed situation and therefore 
the time evolution operator of the system, after the interaction is switched on, writes U(f, t) = 
ex.p(—iVA A (r)t/%) and with this tp(r,t) = U(r,t)i(>o(r). Here we demand that the density 
distribution does not change during the interaction time, which is legitime in the so called 
Raman Nath regime. The Raman-Nath approximation is vàlid as long the gained kinetic 
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energy is much smaller than the interaction potential. The momentum distribution after the 
interaction time t is given by 



n(k, t) 
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To evaluate the magnitude of the dipole moment, we consider the àtoms as two level 
systems driven by a coherent light ficld. This system has a steady state dipole moment, which 
arises from a superposition of the atòmic ground and excited state. The amplitude of this 
oscillating dipole can be derived using optical Bloch equations [17]. The steady state dipole 
moment reads 



d = 2- 9e 
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where d ge is the dipole matrix element, = TyJ Iq/2I ss x the Rabi frequency, T the natural 
linewidth, 8 the làser detuning from resonance and the saturation parameter 
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lo is the intensity of the flashing làser and 7 sa t = irhcT/3\ 3 the saturation intensity. A 
maximum dipole moment only depends on atòmic parameters 
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and is reached for any detuning if the intensity is set to Iq = I sa t(l + 4i5 2 /r 2 ) with ojo being 
the transition frequency. The maximum dipole moment coincides with a saturation parameter 
of one and therefore the spontaneous scattering rate is fixed to T/4. 




Fig. 1 - A Bose-Einstein condensate within an one-dimensional optical lattice. The resulting density 
distribution is a stack of pancake-shaped condensates. A làser along the y-direction induces the 
dipole-dipole interaction. The linear polarization of the flash beam can be altered from x-polarization 
(ip=0°) to z-polarization (ip=90°). The altered momentum distribution is either projected onto the 
x-z plane or the y-z plano. 

Numerical calculations-realistsic experimental setups. - To illustrate the mechanical ef- 
fect of the dipole-dipole interaction, we choose a specific example which corresponds to pa- 
rameters that are realistic in a typical Bose-Einstein condensate experiment with Rb 87 àtoms 
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in the F=2, mp=2 ground state. The atòmic cloud is confined in a cigar-shaped magnètic 
trap with the long axes along z as depicted in fig. ^ In addition, we adiabatically switch on a 
retroreflected làser beam at 785nm to create an optical lattice along the z-axis to increase the 
density. The depth of the lattice is set to 100 recoil energies, which results in an axial trapping 
frequency of 105 kHz. The additional axial confinement due to the magnètic trap is neglected. 
Radially, the trapping frequency is set to 1 kHz. Due to the strong axial confinement, the 
density distribution of the ground state can not be calculated in the Thomas-Fermi approxi- 
mation. We solved the full Gross-Pitaevskii equation [18] numcrically with an imaginary timc 
Schròdinger equation for 250 àtoms in a single lattice site for the given trapping freqüències. 
The resulting pancake-shaped density distribution can be approximated in radial direction by 
a paràbola with a Thomas-Fermi radius of 1.15 /zm and axially by an Gaussian distribution 
with a 1/e 2 radius of 34.2 nm. The resulting peak density is 9.7 x 10 20 m -3 and the chemical 
potential is 5.8 kHz. In the following, we assume an infinite stack of equal pancakes separated 
by A/2=785/2 nm. The calculation is carried out initially on a single pancake and is finally 
superposed according to the infinite stack. The grid for the numerical calculation is chosen 
to be 64 x 64 x 128 = 65.536 lattice points and the grid lattice spacing is A ;/32 = 24.5 nm 
where A; is the wavelength of the làser generating the optical lattice. The direction of the 
dctuning of the làser affects only slightly the lattice constant and with this the superposition 
of the constituent potentials of each pancake. 

The intensity of the flash beam is set to 1120 I sa t and its frequency is detuned 100 MHz 
from the F=2 to F=3 transition (D2 multiplet of Rubidium) at 780.249 nm, which corresponds 
to a detuning of 16.7 T. Using such a large detuning, one can neglect an inhomogeneous 
illumination of the atòmic cloud since only a small fraction of the light is absorbed. The 
induced dipole potentials alters the effective detuning to the flash beam which results in an 
altering phase and magnitude of the oscillating dipoles within the cloud. This effect can also 
be neglected, since the detuning of the flash beam is much larger than the induced dipole 
potential. Also nonlinear effects as lensing by the inhomogeneous density distribution and 
radiation trapping are strongly suppressed. The flash beam propagates along the y-direction 
and its polarization angle <ç can be altered from 0° (polarization along the x-axis) to 90° 
(polarization along the z-axis). The steady state dipole moment with this parameters is 5.26 
Debye. The flash time was set to 300 ns, long compared to 1/r, so that all dipoles are in 
a steady state, but short enough for being in Raman Nath regime. This means also that 
superradiant effects can be neglected [19]. The àtoms are initially prepared in the F=2,mi?=2 
state in respect to the z-axis. 

To account for additional broadening effects of the momentum distribution by spontaneous 
scattering events, the time evolution of the full atòmic density matrix was carried out. This 
includes the different polarizations of the driving field, the pumping of the àtoms into other 
mp states and the angular distribution of the aequired recoil. For given parameters, about 
three photons are scattered per atom. The amount of scattered photons can be experimentally 
checked for consisteney purposes by the shift of the center of mass position after some timc 
of flight of the atom cloud. 

The final momentum distribution is given by the convolution of the momentum gain due 
to the induced dipole potential, the mean field of the Bose-Einstein condensate and the spon- 
taneously scattered photons. To extract a mean momentum broadening, the convoluted dis- 
tribution was fitted with very good agreement by Gaussian distribution. 

Results. - In fig. [3 the distribution of the light-induced dipole-dipole interaction po- 
tential through the center of a pancake along the y-axis is shown. The contribution of the 
neighbouring pancakes is included by periòdic boundary conditions. The different curves rep- 
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Fig. 2 - The left figure shows the induced potentials along the y-axes at x=z=0 for difïerent radial 
widths p of the atòmic clouds. The numbers and the vertical bars at each plot indicate the radial 
Thomas-Fermi radius. The peak atòmic density was kept constant at 9.7 x 10 20 m -3 for all widths. 
The polarization of the flash beam points always along the x-axis. On the right side, the maximum 
acceleration as dependence of the radial Thomas-Fermi radius p is depicted. The data shows the 
maximum gradient of the potentials divided by the mass of a Rubidium atom (m=1.44 xlO~ 25 kg). 



The horizontal line a) gives an upper limit for the unidirectional acceleration {—$-/rn = hkT/Am) 
due to radiation pressure. The horizontal line b) marks the maximum acceleration due to momentum 



diffusion caused by spontaneous emission processes (- 
100 MHz and an intensity of 1120 I sat . 



-/m). The flashing làser has a detuning of 
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Fig. 3 - The two graphs depict the width of the momentum distribution of the atòmic cloud after 
the light-induced dipole-dipole potential was applied for 300ns. The Thomas Fermi radius was set to 
1.15 /im. The graph on the left shows the broadening in momentum space projected onto the y-axes 
as a function of the polarization angle and the right the projection onto the x-axes. The slash-dotted 
line includes the effect of the light-induced potentials and the chemical potential. The sòlid line is the 
incoherent background of the spontaneous scattered photons as a function of the polarization angle 
ip. To receive the full width in momentum space one has to add the two curves in quadraturc. 
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resent pancakes with different radial sizes but at fixed axial size and constant density and show 
the dependence of the induced potentials on the geometry of the atòmic cloud. In the limit of 
an infinite cloud with a constant density distribution, the potential would be just a constant 
and its mechanical effect on the àtoms vanishes. The increase of the induced potentials in the 
center of the cloud with larger radii arises form the greater atom number within the cloud, 
since the density is kept constant. The interaction potential is on the order of several MHz, 
which is large compared to all other energy scales in the system like the trapping freqüències 
and the chemical potential due to the interaction via s-wave scattering. 

The right hand side of fig. [SJshows the maximum acceleration extracted from the potentials. 
For sufficiently small radial sizes of the atòmic cloud, the maximum acceleration prevails the 
unidirectional acceleration due to the spontaneous light forcc. This allows to clearly 
distinguish the effects emerging from the dipole-dipole interaction from spontaneous scattering 
events. 

Finally the slash-dotted curves in fig. [3] show the calculated widths of the momentum 
distribution for the previous parameters along the x and y direction as a function of the 
polarization angle. Not includcd in the slash-dotted curves is the contribution of the sponta- 
neous scattering events represented by the sòlid lines. The broadening in momentum space 
can be up to 10 recoils, which is fairly larger than the contribution of the chemical potential. 
Noticable is the existence of a strong dependence of the broadening on the polarization angle 
ip. In both directions exists an angle at which the effect of the dipole-dipole interaction nearly 
vanishes and the momentum distribution is dominated by the released chemical potential. 
The plot on the left side shows a minimum close to the so-called màgic angle at 54.74° where 
the interaction of two dipoles vanishes. Such a minimum broadening is a clear signature of the 
dipolar character of the potentials since it can not be explained by other light-atom interaction 
mechanisms. 

Conclusion. - We have identified a new regime of coherent light-atom interaction, where 
novel coherent mechanical effects due to dipole-dipole interactions are predicted. The cal- 
culations show that this mechanical effect of the light-induced potentials can be detected 
experimentally for realistic experimental parameters. The distinct angular dependence of the 
polarization of the light field is a clear indication for a dipole-dipole interaction. The numer- 
ical treatment of the incoherent background of spontaneous scattered photons show that the 
signal of the desired mechanical effect is about ten times larger and exhibits not such a char- 
acteristic angular dependence. By carrying out the proposed experiment, one gains insight in 
the physics of coherent dipole-dipole interaction and the feasibility for its usage for quantum 
information processing [10]. 

We thank M. Lewenstein and L. Santos for helpful discussions and the DFG for the financial 
support within the priority program SPP1116. 
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